TO APPEAR IN THE JULY 2002 ISSUE OF THE ASTRONOMICAL JOURNAL 

Preprint typeset using IATf^X style cmulateapj v. 14/09/00 



< 



(N 



(N 
^1" 



Oh 



c3 



SEARCHING FOR PLANETS IN THE HYADES II: SOME IMPLICATIONS OF 
STELLAR MAGNETIC ACTIVITY 1 

Diane B. Paulson 

Astronomy Department, University of Texas, Austin, TX 78712 
apodis@astro.as.utexas.edu 

Steven H. Saar 

Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
ssaar@cfa.harvard.edu 

William D. Cochran 



. McDonald Observatory, University of Texas, Austin, TX 78712 

o . 

O AND 

<N 

, . Artie P. Hatzes 

^ ' Thiiringer Landessternwarte Tautenburg, D-07778 Tautenburg, Germany 

To appear in the July 2002 issue of The Astronomical Journal 



ABSTRACT 

The Hyades constitute a homogeneous sample of stars ideal for investigating the dependence of planet 
formation on the mass of the central star. Due to their youth, Hyades members are much more chro- 
mospherically active than stars traditionally surveyed for planets using high precision radial velocity 
techniques. Therefore, we have conducted a detailed investigation of whether magnetic activity of our 
Hyades target stars will interfere with our ability to make precise radial velocity (i> ra d) searches for 
substellar companions. We measure chromosphcric activity (which we take as a proxy for magnetic 
activity) by computing the equivalent of the i? HK activity index (which is corrected for photospheric 
contributions), from the Ca n K line. (-Rhk) i s n °t constant in the Hyades: we confirm that it decreases 
| with increasing temperature in the F stars, and also find it decreases for stars cooler than mid-K. We 

examine correlations between simultaneously measured i? HK and radial velocities using both a classical 
statistical test and a Bayesian odds ratio test. We find that there is a significant correlation between 
i? HK and the radial velocity in only 5 of the 82 stars in this sample. Thus, simple i? HK ~~ u rad correlations 
will generally not be effective in correcting the measured v ra d values for the effects of magnetic activity 
in the Hyades. We argue that this implies long timescale activity variations (of order a few years; i.e., 
magnetic cycles or growth and decay of plage regions) will not significantly hinder our search for planets 
in the Hyades if the stars are closely monitored for chromosphcric activity. The trends in the radial 



velocity scatter (a' v ) with (-R HK ), vsini, and P TOt for our stars is generally consistent with those found 
in field stars in the Lick planet search data, with the notable exception of a shallower dependence of a' v 



on (-R^k) f° r F stars. 



' Subject headings: clusters: open (Hyades) — stars: planetary systems — techniques: radial velocities 

— stars: chromospheres — stars: activity 

1. introduction neous sample of target stars. By varying only one stel- 

m, ,. , .. , .. / / lar parameter at a time, we may begin to understand the 

The high-precision radial velocity u rad technique e.g., v ,. , , • , , , r 

° , /-mnc^ u u i ui c . i ■ properties and processes which govern planet formation. 
Butler et al. (1996)) has been remarkably successful m TT ^ , , , . r ?,.,,„, . 

. , , , i We have been searching lor planets m the Hyades clus- 

discovermg planetary-mass companions to nearby mam- . , „„„ , . , ,„„„,> , r , n T \ 

? , , . , U1 I u ■ ter since 1996 (Cochran et al. (2001), hereafter Paper I), 

sequence stars. Studying these systems enables us to begin _ . . v , . , v " . . f / 

. , , , . , ,i n , r Ihe Hyades cluster provides a uniform sample ol stars m 

to place important new constraints on tne planet forma- . „. . „, „ . . . . 

c c , ■ , . . j- i i -j. age and initial metallicity. Stellar mass is tne mam mde- 

tion process. So tar, high precision radial velocity surveys b , ... . _. . r . 

1 . . , , ■ i,„„, j v ,i pendent variable among the members, thus, it planetary 
have primarily targeted the brightest, and hence the near- * , r . <3 ... ... .. , , , 

, / ; . . ,. , , , , • , , systems are found, we will easily be able to study the role 

est cnromospnerically inactive, slowly-rotating solar-type ,. , • , , r , • 

i. rru. v.-i ii, i a -a c stellar mass plays m the planet formation process, 

stars. Ihus, while the targets cover a wide range ol ages, T ,,,,,,, , 

... j , , . , . ,v j ■ , i In order to detect planetary systems, the w ra d signa- 

compositions, and perhaps histories, they are dominated , . . . . J , J 

, , . , , ture must be unambiguous, and we must be certain it 

by old, higher mass stars. . , , . ., . p . ' . 

A i ±_ j ; j- t j-i, i j- t i- ' is due to an orbiting body and not intrinsic to the star. 

A more complete understanding ol tne planet formation _ . . ° . J , ,. r . .. 

u . , , . . , . , , ,. Several groups have explored the implications ol intrinsic 

process lunges upon determining which stellar properties , ., ° . . , . , „ , . m „ 

■ a , , , * ,. j , . Tl . stellar properties, m particular stellar pulsations (Gray & 

influence planetary system formation and dynamics. Ihis \ j 

can be accomplished with a well controlled and homoge- 

1 Data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California 
Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made 
possible by the generous financial support of the W.M. Keck Foundation. 
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Hatzcs 1997), rotational velocity, and chromospheric activ- 
ity (Saar & Donahue 1997; Saar et al. 1998; Santos et al. 
2000a) on radial velocity measurements. Stellar non-radial 
(as well as radial) pulsation can easily give observable ra- 
dial velocity variation, but may be distinguished from or- 
bital motion by their resulting variations in stellar line 
bisectors (Gray & Hatzes 1997; Hatzes et al. 1997; Hatzes 
et al. 1998a, b; Brown et al. 1998b, a). The rotational mod- 
ulation and evolution of velocity structures or non-uniform 
active regions can cause similar u ra d signatures to orbiting 
planets. Walker et al. (1992) showed that in the case of 7 
Cep, the w ra d measurements suggestive of a planet were in- 
stead likely due to rotation of active regions seen in equiv- 
alent widths of the Ca 11 8662 A line. Walker et al. (1995) 
detected a correlation between w ra( j and long-term Ca II 
emission variations in k 1 Ceti probably due to its mag- 
netic cycle. Queloz et al. (2001) showed that the apparent 
radial velocity signal in HD 166435, which was originally 
interpreted as due to a planetary companion is actually 
correlated with photospheric line profile variations and is 
best explained by the presence of dark photospheric spots. 
Dravins & Nordlund (1990) discussed the possibility of 
apparent false Doppler shifts over the course of a stellar 
activity cycle, though, McMillan et al. (1993) has shown 
that solar cycle magnetic variations are undetectable in 
disk integrated spectra of the sun. Saar and collaborators 
(Saar et al. 2002; Pearson et al. 2000; Saar & Fischer 2000; 
Saar & Snyder 1999) are working on methods for correct- 
ing stellar radial velocities to compensate for the stellar 
activity, but these techniques are still under development. 
Because of worries about such activity-related v ia d noise, 
past planet surveys have been biased against young, active 
stars (Vogt et al. 2002; Henry et al. 1997; Saar & Donahue 
1997). Consequently, out of 76 planets discovered as of 
December 2001, only four have been found around young, 
active stars: e Eridani (Hatzes et al. 2000), l Horologii 
(Kiirster et al. 2000), HD 192263 (Santos et al. 2000a), 
and GJ 3021 (Naef et al. 2001). 

Because the Hyades cluster is so young (625 ± 50 Myr; 
Pcrryman et al. (1998)), the stars can be quite active. 
Thus, we must determine if this activity is a major con- 
tributor to the observed w ra d variances. In this paper, we 
explore the connection between activity variations and ra- 
dial velocity measurements in Hyades stars by looking for 
correlations between the Ca 11 K emission and the v ra d 
(measured in the same spectra). This will help us un- 
derstand to what degree this phenomenon masks or com- 
plicates detection of the low velocity signal induced by 
planetary mass companions. 

2. OBSERVATIONS AND ANALYSIS 

We are conducting a radial velocity survey of Hyades 
dwarfs using the Keck HIRES spectrograph to search for 
Jupiter and Saturn-mass companions (see Paper I). By us- 
ing a 10m class telescope, we are able to study even low 
mass Hyads (down to M m 0.4M Q ). So far, this survey has 
produced a number of stars showing rms velocity variation 
significantly larger than the internal velocity uncertainty 
(<7j, which includes instrumental noise, photon noise, and 
increased noise due to observational effects such as reduced 
line density and depth in F stars and increased line width 
in stars with higher v sin i) . Unfortunately, the scheduling 



of our observing time on Keck has made it impossible to 
sample short periods without significant aliasing. We be- 
gan observations of short period candidates with the High 
Resolution Spectrograph on the Hobby-Eberly Telescope 
in the fall of 2001. 

2.1. Sample 

Our total sample consists of 98 Hyades dwarfs ranging 
from F5 to M2. All of the targets are confirmed mem- 
bers according to Perryman et al. (1998), and all known 
spectroscopic binaries have been removed. We imposed a 
stellar rotational velocity limit of vsini < 15 km s _1 . For 
the chromospheric activity study presented here, we ne- 
glect all stars with fewer than six observations, reducing 
the sample to 82 stars. 

2.2. w ra d measurements 

The observations at Keck I make use of the HIRES (Vogt 
et al. 1994) with I 2 gas absorption cell as a standard ve- 
locity reference (Valenti et al. 1995) as detailed in Paper 
I. S/N~300 is required in this configuration (at resolution 
R ~ 60,000) to obtain v ra d precision (in the best cases) of 
^3 m s _1 (see Paper I); S/N^100 yields precision of ~5-6 
m s _1 . In addition, the exposure times must be kept to 
<15 minutes. We use standard IRAF packages to reduce 
the CCD images and extract the observed spectra. 

It is useful to explore how the rotational broadening in 
the spectra affects the internal v ra d errors (<7j). Since we 
have preselected stars with vsini < 15 km s _1 , we only 
have a small range of v sini to study. In Figure 1, we 
show the relationship between v sin i (taken from the lit- 
erature as referenced in Table 1) and the mean internal 
errors ((oi)) we achieve. Stars with only vsini upper lim- 
its have not been included in the analysis. Error bars on 
(<7j) indicate the rms about the mean. We note that even 
with this limited set of data, we see the trend that one 
would expect; that increasing vsini degrades our velocity 
precision, due to the lowered precision in determining the 
center of broader lines. In order to estimate how large our 
errors are at vsini=15 km s _1 , we fit a linear relation to 
the data. We find 

(crj) = 1.050 + 0.555* v sin i m s" 1 . (1) 

Therefore, at the cutoff of 15 km s _1 , we can still achieve 
Vrad precision of <~ 10 m s _1 . We would like to point out 
that these are our mean internal errors of all of our spec- 
tra of a given star. If instead we had chosen the highest 
S/N spectrum, we would have derived Ui = 6.5 m s _1 at 
vsini = 15 km s _1 . For vsini < 3 km s _1 , other sources 
of line broadening become important (e.g., macroturbu- 
lence), and (<7j) is no longer dominated by vsini. 

2.3. Calibration of Chromospheric Activity 

Most stellar chromospheric activity is attributed to the 
interaction of magnetic fields with convection (Middelkoop 
& Zwaan 1981; Tinbergen & Zwaan 1981; Middelkoop 
1982). However, the effect of magnetic processes on the 
integrated radial velocity of a star at the few m s -1 level 
is only beginning to be explored. To date, stars which are 
chromosphcrically active and show no obvious trend in the 
velocity signal over a given timescale are either observed 
less frequently or removed from planet surveys (Vogt et al. 
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2002; Cumming et al. 1999; Saar & Donahue 1997). This is 
because these stars tend to show somewhat higher levels of 
radial velocity scatter, attributed to the magnetic activity, 
and this will possibly inhibit detection of low-amplitude 
velocity signals. 

The question arises as to whether stellar activity in 
Hyades stars induce significant periodic centroid shifts in 
photospheric absorption lines, and thus v rac i signals which 
could be confused for perturbations made by planetary 
companions. The members of the Hyades cluster are chro- 
mospherically active, and quite young. Because activity 
might be a significant source of radial velocity variation 
(Saar & Donahue 1997; Saar et al. 1998; Saar & Fischer 
2000), we include the Ca n H & K lines in the spectral 
region of each exposure for velocity measurement with the 
Keck HIRES. We have monitored the chromospheric ac- 
tivity of our target stars, acquiring Ca n K emission core 
flux measurements simultaneously with each velocity mea- 
surement. 

To measure stellar chromospheric activity, the Mt. Wil- 
son S index is adopted. This index is defined (e.g. Bali- 
unas et al. (1995)) as a quantity proportional to the sum 
of the flux in lA FWHM triangular bandpasses centered 
on the Ca II H and K lines divided by the sum of the flux 
in 20A bandpasses in the continuum at 3901 and 400lA 
(Soderblom et al. 1991; Baliunas et al. 1995). 

At Keck 1, we set HIRES so that the Ca II H and K lines 
are contained within the spectra for all stars observed. The 
four quantities to be measured (the two calcium line core 
fluxes plus the two continuum bandpass fluxes) are spread 
across three echelle spectral orders which overlap by a few 
A. As a result of flat fielding uncertainties in the most 
blucward order, we do not use measurements in this or- 
der, i.e. of the blue (390lA) continuum. In addition, we 
do not measure the Ca n H line (at 3968. 47A) because 
for some stars, the wings of strong Balmer He features (at 
3970. 07A) lay within the measured bandpass. This would 
affect measurement of the Ca II H line flux as the absorp- 
tion can be strong enough so as to artificially decrease the 
Ca II H flux, and strong He emission (seen in a few of the 
M dwarfs) will artificially increase the measured Ca n H 
line flux. Therefore, we have defined an index S'Keck which 
is the ratio of the flux in a lA triangular bandpass centered 
on the Ca n K line to the flux in a 20A bandpass centered 
on the redward continuum at 4001 A. Figure 2 shows the 
numerical filter we used in measuring the Ca n K flux. 
Several of our Hyades program stars have previously been 
measured as part of the Mt. Wilson program (Duncan 
et al. 1984). In Figure 3, we show the relationship be- 
tween our -S'Keck and Sbuncan (the published mean S data 
for stars measured by Duncan et al. (1984) including error 
bars which characterize the rms of the measured S index 
over the course of observations (including variability). We 
find the linear relation 

^Duncan = -0.027(±0.046) + 0.991 (±0. 120) S K cck (2) 

Both the slope and the intercept of the formal fit are within 
1 a of S'Keck = SWncan- We are * nus t° transform 
our data into a standard Mt. Wilson S index scale using 
Equation 2. We call our transformed values Smw- 

In addition, we further transform our chromospheric S 
index to i?^ K = F^ K /F bo i, where is the Ca n HK 
surface flux, corrected for flux contributions from the pho- 



tosphere, and -Fboi is the bolometric flux (Noyes et al. 1984; 
Middelkoop 1982). Thus, R' HK is the normalized, purely 
chromospheric component of the HK flux and is also useful 
as a means to compare stars of different spectral types. To 
convert from Smw to -Rhk, we used the method outlined 
in Noyes et al. (1984), with B — V color index of each star 
obtained from Allende Prieto & Lambert (1999), and no 
reddening corrections were used. 

3. RESULTS 

3.1. Chromospheric activity 

Figure 4 shows the variation of (-R^k) w ith B — V for 
our sample of Hyades dwarfs. This demonstrates that ac- 
tivity as measured by (i?n K ) i s n °t constant with B — V in 
Hyades stars: (-Rhk) values decrease from a plateau with 
(logi^K> w -4.42 for B - V <0.7 and for B - V >1.1. 
Noyes et al. (1984) notes that the values of R' HK for 
B — V >1.0 may be uncertain by as much as 20%. How- 
ever, the downward trend of B — V below 0.7 has not 
yet been addressed to our knowledge. A few outliers with 
B — V <0.7 (see Fig. 4) were determined not to be Hyades 
members by D. Latham (private communication) based 
on too large Hipparcos distances, photometry below main 
sequence and velocities which disagree with members, in 
contrast to Perryman et al. (1998) who found these same 
stars to be members. The decline for B — V <0.7 confirms 
the trend observed by Duncan et al. (1984) and Soderblom 
(1985) and seen in the data collected by Rutten (1987). 
HK data for the younger Pleiades cluster (age ~ 0.1 Gyr) 
shows more scatter (Rutten 1987), but a similar, though 
weaker, decline for B — V <0.7. This suggests that the 
decline in activity for B — V <0.7 is real (as opposed to 
a problem in R[i K (B — V) calibration); a color-dependent 
calibration error should not evolve significantly between 
the ages of the Pleiades and Hyades. Thus, although we 
must be cautious about stars with B — V >1.0, we confirm 
that at fixed age, there is a broad maximum in chromo- 
spheric emission for 0.7< B - V < 1.0 (wG5 to K3), with 
decreases for both hotter and (tentatively) for cooler stars. 

3.2. Intrinsic u ra( j Jitter 

3.2.1. Comparison of Observed v rac i Jitter with Empirical 

Relations 

Saar et al. (1998) (hereafter SBM) used the observed 
radial velocity dispersion (o~ v ) in the Lick planet search 
data to compute the v rac i dispersion above the internal 
errors, a' v — (of — of) 05 . They then searched for cor- 
relations between a' v and (-Rhk)' rotational period (P TO t), 
and projected rotational velocity (v sini). They also devel- 
oped a simple model based on rotating/evolving starspots 
and plage to predict a' v . The Hyades dwarfs studied here 
are significantly more active on average than those stars 
studied in SBM, making it useful to investigate whether 
Hyad o~' v behave similarly. We calculated o~' v in the same 
manner as SBM, removing internal noise from each mea- 
surement. Due to the selection of stars with vsini <15 
km s _1 (Paper I), we have also selected for stars with low 
radial velocity jitter as shown in the model developed in 
SBM. We have also ignored stars with very large (3> 100 
m s _1 ) systematic velocity trends. F stars in our sample 
(ignoring one outlier) have a' v ranging from ^8 to 38 m 



4 



s -1 with an average value of (cr' v ) — 14.8 ± 8.4 m s -1 , G 
stars range from ~5 to 40 m s -1 ((<r' v ) = 18.9±9.0 m s^ 1 ), 
K stars (ignoring one outlier) range from ^3 to 45 m s _1 
{(a' v ) = 14.2 ± 13.6 m s" 1 , and M stars from -4 to 20 m 
s_1 (K) = 12.7 ±7.6 m s- 1 ). 

A comparison of the observed a' v and the values pre- 
dicted by empirical relationships found in the Lick v rac i 
database (cr^ prcd ' ) ; see SBM) are shown in Fig. 5 (first 
three panels) and summarized in Table 2. For the rotation- 
related quantities, P Iot and vsini, agreement is generally 
quite good, with the scatter a about the line o~' v = cr^ prcd ' ) 
at or below the values found for the empirical SBM fits. 
(We have considered only stars with measured P ro t ) Stars 
with only upper limits to vsini are also all in agreement 
with predictions. 

Only relatively few Hyades stars have measurements 
of P ro t and vsini needed to compare a' v with the pre- 
dictions of cr^ prcd - ) based on these parameters using the 
SBM results. With (i?n K ) though, all the Hyades sam- 
ple can be used except M stars (which were not stud- 
ied in SBM). If cr^ prcd ' is predicted using i?^ K and the 
combined Lick G and K star sample, the scatter a rms of 
Hyades G and K stars' a' v about this o$ pred ^ (i?^ K (G+K)) 
is slightly higher than for the Lick stars. This is true even 
if one Hyades outlier (a K star) is removed (a rms (Hyades) 
= 0.28 dex compared with <r rrns (Lick) = 0.25). To ex- 
plore this result further, we compared a' v for the Hyad 

G and K stars separately with cr^ prcd - ) (i?^ K (G+K)), and 

also with av Wcd ^ based on separate empirical fits to the 
Lick G and K stars (not published in SBM). These tests 
indicate that a rms (Hyades) is smaller than ay TOS (Lick) 
for the G stars («0.22 dex), whether a' v is compared 

to CT ; (prod) (i^ K (G+K)) or a; (prod) (i^ K (G)). The scat- 
ter is notably worse for K stars, though improved some- 
what with the separate K star empirical fit: o~ rms «0.32 
dex using o-u Pred ^(i?H K (G+K)) , o rms w0.28 dex using 
^./(pred) (iJ^ K (K)), with one outlier removed in each case. 
Thus, since F, G, and K stars are best fit when treated 
separately, the Hyades a' v data suggest a systematic spec- 
tral type dependence in the relationship between a' v and 

The Hyades F stars (excluding one outlier) show a sys- 
tematic trend relative to the a' v = cr^ prcd ' ) line which sug- 
gests the empirical relation (cr' v (F)oc (R^) 1 ' 7 ; SBM) may 
be too steep. Indeed, the scatter is lower and the trend 
largely removed if the relations for G or G+K stars are 
used instead (a rms = 0.21 dex or 0.28 dex, respectively, 
compared with cr rms = 0.38 dex using the F star fit). The 
vsini < 15 km s _1 limit of our sample, though, may be a 
factor; while this limit has little effect on G and K Hyades 
stars, it definitely excludes some high vsini (and hence 
high R' HK ) F stars. The Lick analysis was performed with- 
out any vsini restrictions, and includes two F stars with 
vsini > 15 km s _1 . Restrictions on vsini will tend to 
limit o' v for a fixed i?HK i smce the higher v sin i is an im- 
portant component for enhanced v ra d noise (Saar & Don- 
ahue 1997). 

The two "outlier" stars noted above (one F and one 
K type) show a' v values significantly enhanced (by > 



2oy ms (Lick)) over that predicted by o~' v (-Rhk)- Hence, 
these stars have large additional v ra d variations not asso- 
ciated with activity, and are particularly good candidates 
for further and more frequent observations in search of low 
mass companions. 

3.2.2. Comparison of Observed v rac i Jitter with a Simple 

Model 

We have also compared Hyades a' v values with the com- 
bined spot/plage a' v model outlined in SBM. This model 
combines the rotating spot model from Saar & Donahue 
(1997) with a (very) simple plage model which assumes 
that convective velocity changes due to plage magnetic 
fields arc cx v mac , the macroturbulent velocity The full 
model is given by 

< |mod) ~ A /[4.6/0- 9 («sini)co S <e)]2 + [a J 4(/ f >)(D mac + Dsini)]2[ms- 1 ]. 

(3) 

Here, f$ is the differential spot filling factor (the portion 
of the total spot filling factor responsible for photometric 
variations). Also, (8) is the mean spot latitude, a is an 
adjustable constant, and A(fp) is a function of the plage 
filling factor fp, given by A(fp) = 0.25 sin(ir ^/Jp) (the 
scaling factor was not explicitly given in SBM; it scales 
the maximum of A(fp) = fp). Following SBM, we take 
f P = 0.08(r c /P rot ) 1 - 8 < 0.65 (Saar (1996); where r c is 
the convective turnover time from (Noyes et al. 1984)), 
(9) = 45°, a s=s 9, and estimate v mac from relations for 
active stars in Saar & Osten (1997). 

Unfortunately, only five stars had the required data 
(vsini, -P ro t, and photometry sufficient to estimate fs) 
to make a complete model estimate of a' v . If we assume 
fgKiO for F stars (Hyades F stars show little photometric 
variability - Radick et al. (1998)) used them as constant 
check stars), and estimate their P ro t from R' HK (vis. Noyes 
et al. (1984)) we gain two more stars. Several more have 
upper limits to vsini. As shown in Figure 5 (lower right) 

and Table 2, the agreement between o~' v and cr^ mod -' for the 
(few!) stars we can model adequately is quite good: o~ rms 
is ~6 to 7 m s _1 or w 0.13 dex (comparable to SBM). 
All cases where only upper limits on a' v (mod) are possible 
(due to v sin i upper limits) are in agreement as well. 

3.3. R' HK vs. v rad 

We have calculated the linear correlation coefficient (r) 
in the standard way (Bevington & Robinson 1992) to 
search for any statistical correlations between the chro- 
mospheric activity and the radial velocity for each star. 
The values of r (which range from -1 to 1) calculated are 
shown in column 6 of Table 3. We then use this quantity 
to test statistically (both in a frequentist and a Bayesian 
manner) whether any linear correlations exist between the 
v rad and i?HK- 

We first follow the frequentist approach given in Bev- 
ington & Robinson (1992) to determine a probability that 
the Vrad and R' HK are correlated. The correlation coeffi- 
cient is used to find a probability P c that the data come 
from an uncorrelated parent population (P c = l indicates 
that the data are completely uncorrelated, while P c =0 in- 
dicates be a completely correlated data set). The results 
from this test are in column 8 of Table 3. In using this 
method, we find several stars with high probabilities that 
the data come from a correlated parent sample. However, 
this method introduces a bias, since P c is calculated by 
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only considering r relative to an uncorrelated parent sam- 
ple. What is perhaps a better test is to compare the data 
to both correlated and uncorrelated parent samples. This 
way, we can test both hypotheses and make a better de- 
termination of possible correlation in the parent sample. 
To do this, we use a Bayesian approach. 

Our Bayesian analysis uses an odds ratio (K) defined 
by Jeffreys (1961). This incorporates prior knowledge of 
a probability distribution of r of the parent sample. Here, 
we choose a prior distribution to be constant and centered 
on (0,0). We calculate a ratio of the probability that v rac i 
and R' HK are uncorrelated versus the probability that they 
are correlated (i.e., K <Cl indicates a strong correlation, 
K 3>1 the lack of one). For example, for K = 2, there 
are 2:1 odds that the sample is uncorrelated (i.e., a 67% 
chance that the sample is uncorrelated) . Column 9 of Ta- 
ble 3 gives the odds ratios. However, Schmitt (1969) notes 
that the strength of correlations worked out in this man- 
ner is somewhat overestimated: correlation probabilities 
<75% are highly suspicious. In contrast to the frequen- 
tist correlation analysis, with the Bayesian techniques, we 
show that 15 stars have greater than 50% chance of being 
correlated. Of these, 5 show a strong chance of correlation- 
greater than 70%. Due to the small number of observa- 
tions, Bayesian statistics will yield more reliable results, so 
these values are the ones the authors favor. Of the remain- 
ing stars, 40 show slight non-correlation between v ra d and 
R' HK and 11 show slight correlation. So, the overwhelming 
majority of stars are slightly to strongly uncorrelated. Fig- 
ures 6 and 7 show examples of stars with uncorrelated and 
correlated trends, respectively, from the previous analysis. 

4. DISCUSSION 

Our fundamental result is that very few Hyads (5 of 82) 
show significant correlations between simultaneous w rac i 
and i?HK measurements. It is therefore important to re- 
view how stellar activity can alter observed radial veloci- 
ties, and what our result then implies for planet detection 
in the Hyades. In the Sun, Ca n emission primarily reflects 
the surface coverage of plage and active network, since it 
is relatively weak over sunspots themselves (e.g. Linsky 
& Avrett (1970) and references therein), and spot area 
is typically small relative to plage area. Thus, v ra d fluc- 
tuations due to spots (e.g., Saar & Donahue (1997)) will 
not have corresponding AR / HKl and so spots should not 
be a significant contributor to w ra( j - -Rhk correlations. 
Plage/network can generate w ra( j fluctuations in two main 
ways. First, plage is slightly brighter in the continuum 
(by a few %) than the quiet Sun. Assuming this is also 
true for stars, rotation of inhomogencous patches of plage 
will cause traveling enhancements at the local intensity- 
weighted rotational velocity, which will translate into ap- 
parent w ra( j changes on P rot timescales. This effect is com- 
pletely analogous to the one caused by spots. Since the 
brightness enhancement of plage is tiny compared with 
the ~90% light deficit (in V) due to sunspots, though, the 
brightness effect of plage on v ra d should be small, even 
considering the typically larger plage area. 

To explore the effects on w rac j-activity correlations, con- 
sider that an identical plage, observed at two rotation an- 
gles ±<f> (measured from disk center), will exhibit identi- 
cal activity enhancements (i.e., AR' HK ((f)) = AR^ K (— <f>)), 



but will show perturbations to v ra d due to brightness of 
opposite sign: Aw rad (0) = -Av rad (-0) = av Iot (<p) (where 
a is some function of the plage brightness enhancement) . 
Thus, a given plage will show a scatter of Av ra d due to 
its brightness as it rotates across the disk. When averaged 
over enough observations, intensity changes due to activity 
(spots or plage) have no net w ra d effect, since the pertur- 
bations oc ±v rot {(f)) will average to zero over the many ob- 
served <fi. Changing the mean level of activity will change 
the average magnitude of the effect and its rms, but once 
again, for sufficient <fi coverage the net Aw ra d = 0. 

In addition to its brightness perturbation, plage also in- 
duces changes in the local velocity field, suppressing con- 
vective velocities in strong magnetic fields. The altered 
velocity field induces changes in the line shape (Livingston 
1982) and leads to an overall convective blueshift in the 
Sun (Cavallini ct al. 1985). Furthermore, these changes 
will vary from line to line based on their strength and 
excitation, just as line bisectors do (e.g., Asplund et al. 
(2000)). Thus, in the case of plage the w ra d perturbation 
may be written Av rad (</>) = v U u C (4>) + v shapc (cf>), where 
Wbiuo(</ ) ) is the overall convective shift of the line core, and 
v s hape(</>) is the apparent w ra d change due to the altered 
line profile shape. The latter arises because most methods 
measure the t> ra d of individual exposures by comparison 
to a single high S/N "template" spectrum (either an av- 
erage spectrum or a single deep exposure). Fluctuations 
in line shape relative to this "template" will inevitably 
lead to an apparent shift in the line centroid, and hence 
w ra d- In contrast with brightness perturbations, the time 
average of w b iuo(0) + Wshape( ( / ) ) will in general be non-zero, 
since both, due to their intimate connection with convec- 
tion are (unlike f ro t) symmetric about cf> — in the time- 
averaged sense. Clearly, strong correlations between v ra d 
and R' HK will result primarily from these long timescale 
changes in average plage area (see also Saar & Fischer 
(2000)). The present results say little about short term 
(timescales < P ro t) changes in w ra d, or about the effects of 
starspots. 

The lack of many significant w ra d - -Rhk correlations in 
the Hyades thus implies that longterm changes in plage ac- 
tivity have little effect on w ra d for our stars. This is not an 
entirely surprising result. If we look at active stars in Bal- 
iunas et al. (1995), relatively few show clear cycles in Ca 
II (see also Saar (2002)). Without a systematic long-term 
activity variation, active stars are less likely to show strong 
w ra d - R'ylk correlations. Typically, short term (rotational) 
jitter in v ra d, and perhaps flares in activity, are expected 
to dominate active stars without strong cycles, swamping 
potential v ra d - R' aK correlations with rapid fluctuations 
in one or both variables. Thus, the method suggested by 
Saar & Fischer (2000) for correcting w ra d timeseries for 
some of the jitter induced by activity will not be effective 
for most Hyads. Other methods, based for example on 
variable line bisector changes (Saar et al. 2001; Saar et al. 
2002; Queloz et al. 2001) may be useful in diagnosing (and 
possibly correcting) for short term u rad jitter from spots 
and plages. 

Clearly, the changing (-Rhk) with B — V implies that the 
inverse Rossby number Ro^ 1 = Tc/P ro t (at least as de- 
fined by Noyes et al. (1984)) is not constant with mass at 
fixed age. This confirms and extends (to B — V > 1.0) the 
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similar conclusion of Soderblom (1985). Consistent with 
this, the age calibration of Donahue, Dobson & Baliunas 
(1997), which estimates age t as a function of (-R H k) a l° ne , 
predicts (logi[yr]) = 8.89 ± 0.27 (in agreement with Per- 
ryman et al. (1998)) but overestimates t for B — V < 0.60 
and B — V > 1.30. The physical implications of a mass 
dependent (-R H k) m the Hyades are less clear. Apparently, 
either generation of magnetic flux, or the physics of chro- 
mospheric heating, or both, vary with mass at fixed age. 
There is some independent evidence for both of these ideas. 
Saar (2001) finds a monotonic relation between magnetic 
flux and i?o _1 , which implies mass dependent magnetic 
flux in the Hyades, given that i?o _1 is not constant. Sup- 
porting the concept of mass-dependent changes in heating, 
several researchers (e.g., Rutten et al. (1989)) have noted 
that in comparison to hotter stars, M dwarfs show Balmer 
emission dominating over HK emission, and chromospheric 
emission in general reduced relative to coronal emission. 

Taken as a whole, the a' v values for the Hyades are in rea- 
sonably good agreement with empirical results and models 



in SBM. The main exceptions are that the Hyades data 
suggest a more continuous change in the dependence of 
the a' v - (-R^k) relationship on spectral type (e.g., G and 
K stars are better considered separately). Also, the SBM 
fit for F stars, a' v oc (-Rhk) 1 ' 7 ma Y be too steep, as also 
suggested by Santos et al. (2000b). The difference may 
however be the result of the limitation of v sinz < 15 m 
s _1 in the present sample. Further analysis of the Lick and 
Hyades data, with and without vsini limits, is needed to 
resolve this issue. 
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and WDC arc supported by NASA grants NAG5-4384 and 
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ported by NASA grant NAG5- 10630. 



REFERENCES 



Allende Prieto, C. & Lambert, D. L. 1999, A&A, 352, 555 
Asplund, M., Nordlund, A., Trampedach, R., Allcnde-Prieto, C, & 

Stein, R. 2000 A&A, 359, 729 
Baliunas, S. L., Donahue, R. A., Soon, W. H., Home, J. H., Frazer, 
J., Woodard-Eklund, L., Bradford, M., Rao, L. M., Wilson, O. C, 
Zhang, Q., Bennett, W., Briggs, J., Carroll, S. M., Duncan, D. K., 
Figueroa, D., Lanning, H. H., Misch, T., Mueller, J., Noyes, R. W., 
Poppe, D., Porter, A. C, Robinson, C. R., Russell, J., Shelton, 
J. C, Soyumcr, T., Vaughan, A. H., & Whitney, J. H. 1995, ApJ, 
438, 269 

Bcvington, P. R. & Robinson, D. K. 1992, Data Reduction and Error 
Analysis for the Physical Sciences (New York: McGraw-Hill) 

Brown, T. M., Kotak, R., Horner, S. D., Kcnnclly, E. J., Korzennik, 
S., Nisenson, P., & Noyes, R. W. 1998a, ApJS, 117, 563 

— . 1998b, ApJ, 494, L85 

Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C, Dosanjh, 

P., & Vogt, S. S. 1996, PASP, 108, 500 
Cavallini, F., Ceppatelli, G., & Righini, A. 1985, A&A, 143, 116 
Cochran, W. D., Hatzes, A. P., & Paulson, D. B. 2001, AJ, submitted 
Cumming, A., Marcy, G. W., & Butler, R. P. 1999, ApJ, 526, 890 
Donahue, R. A., Dobson, A. K., & Baliunas, S. L., Solar Phys. 171, 

211 

Dravins, D. & Nordlund, A. 1990, A&A, 228, 203 

Duncan, D. K., Frazer, J., Lanning, H. H., Baliunas, S. L., Noyes, 

R. W., & Vaughan, A. H. 1984, PASP, 96, 707 
Fckcl, F. C. 1997, PASP, 109, 514 
Gray, D. F. & Hatzes, A. P. 1997, ApJ, 490, 412 
Hatzes, A. P., Cochran, W. D., & Bakker, E. J. 1998a, Nature, 391, 

154 

— . 1998b, ApJ, 508, 380 

Hatzes, A. P., Cochran, W. D., & Johns-Krull, C. M. 1997, ApJ, 478, 
374 

Hatzes, A. P., Cochran, W. D., McArthur, B., Baliunas, S. L., 
Walker, G. A. H., Campbell, B., Irwin, A. W., Yang, S., Kiirster, 
M., Endl, M., Els, S., Butler, R. P., & Marcy, G. W. 2000, ApJ, 
544, L145 

Henry, G. W., Baliunas, S. L., Donahue, R. A., Soon, W. H., & Saar, 

S. H. 1997, ApJ, 474, 503 
Jeffreys, H., S. 1961, Theory of Probability (Oxford: Clarendon 

Press) 

Kraft, R. P. 1965, ApJ, 142, 681 

Kiirster, M., Endl, M., Els, S., Hatzes, A. P., Cochran, W. D., 

Dobcreiner, S., & Dennerl, K. 2000, A&A, 353, L33 
Linsky, J. L. & Avrett, E. H. 1970, PASP, 82, 169 
Livingston, W. 1982, Nature, 287, 208 

McMillan, R. S., Moore, T. L., Perry, M. L., & Smith, P. H. 1993, 

ApJ, 403, 801 
Middelkoop, F. 1982, A&A, 113, 1 
Middelkoop, F. & Zwaan, C. 1981, A&A, 101, 26 
Naef, D., Mayor, M., Pepe, F., Queloz, D., Santos, N. C, Udry, S., 

& Burnet, M. 2001, A&A, 375, 205 
Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & 

Vaughan, A. H. 1984, ApJ, 279, 763 



Pearson, N. S., Saar, S. H., Hatzes, A. P., & Paulsen, D. 2000, BAAS, 
Vol. 197, 1107 

Perryman, M. A. C, Brown, A. G. A., Lebreton, Y., Gomez, A., 

Turon, C., Cayrel de Strobcl, G., Mermilliod, J. C, Robichon, N., 

Kovalevsky, J., & Crifo, F. 1998, A&A, 331, 81 
Queloz, D., Henry, G. W., Sivan, J. P., Baliunas, S. L., Beuzit, J. L., 

Donahue, R. A., Mayor, M., Naef, D., Perrier, C, & Udry, S. 2001, 

A&A, 379, 279 

Radick, R. R., Lockwood, G. W., Skiff, B. A., & Baliunas, S. L. 1998, 

ApJS, 118, 239 
Rutten, R. G. M. 1987, Ph.D. Thesis, 129 

Rutten, R. G. M., Zwaan, C, Schrijver, C. J., Duncan, D. K, & 

Mewe, R. 1989, A&A, 219, 239 
Saar, S. 2001, in Cool Stars, Stellar Systems, and the Sun 11, ASP 

Conference Series Vol. 223, ed. R.J. Garcia Lopez, R. Rebolo, & 

M.R. Zapatro Osrio, p. 292 
Saar, S. & Ostcn, R. 1997, MNRAS, 284, 284 

Saar, S. H. 1996, in Stellar Surface Structure, IAU Symp. 176, ed. 

J. L. K. Strassmeier (Kluwer), 237 
Saar, S. H. 2002, in Stellar Coronae in the Chandra and XMM- 

Newton Era, ASP Conf. Series, in press, ed. F. Favata & J. Drake 
Saar, S. H., Butler, R. P., & Marcy, G. W. 1998, ApJ, 498, L153, 

=SBM 

Saar, S. H. & Donahue, R. A. 1997, ApJ, 485, 319 
Saar, S. H. & Fischer, D. 2000, ApJ, 534, L105 

Saar, S. H., Fischer, D., Snyder, N., & Smolec, R. 2001, in Cool 
Stars, Stellar Systems and the Sun 11, ASP Conf. Series, Vol. 223, 
CD-1051 

Saar, S. H., Hatzes, A. P., Cochran, W. D., & Paulson, D. B. 2002, 
in Cool Stars, Stellar Systems and the Sun 12, ASP Conf. Series, 
in press 

Saar, S. H. & Snyder, N. 1999, BAAS, Vol. 194, 10502 

Santos, N. C, Mayor, M., Naef, D., Pepe, F., Queloz, D., Udry, S., 

& Blecha, A. 2000a, A&A, 361, 265 
Santos, N. C, Mayor, M., Naef, D., Pepe, F., Queloz, D., Udry, S., 

Burnet, M., & Revaz, Y. 2000b, A&A, 356, 599 
Schmitt, S. A. 1969, Measuring Uncertainty: an Elementary 

Introduction to Bayesian Statistics (Reading, Mass.: Addison- 

Wesley) 

Soderblom, D. 1985, AJ, 90, 2103 
Soderblom, D. R. 1982, ApJ, 263, 239 

Soderblom, D. R., Duncan, D. K., & Johnson, D. R. H. 1991, ApJ, 
375, 722 

Strassmeier, K. G., Washuettl, A., Granzer, T., Scheck, M., & Weber, 
M. 2000, A&AS, 142, 275 

Tinbergen, J. & Zwaan, C. 1981, A&A, 101, 223 

Valenti, J. A., Butler, R. P., & Marcy, G. W. 1995, PASP, 107, 966 

Vogt, S. S., Allen, S. L., Bigclow, B. C, Bresee, L., Brown, B., 
Cantrall, T., Conrad, A., Couture, M., Delaney, C, Epps, H. W., 
Hilyard, D., Hilyard, D. F., Horn, E., Jern, N., Kanto, £>., Kcane, 
M. J., Kibrick, R. I., Lewis, J. W., Osborne, J., Pardeilhan, G. H., 
Pfistcr, T., Ricketts, T., Robinson, L. B., Stover, R. J., Tucker, 
D., Ward, J., & Wei, M. Z. 1994, Proc. Soc. Photo-opt. Inst. Eng., 
2198, 362 



7 



Vogt, S. S., Butler, R. P., Marcy, G. W., Fischer, D. A., Pourbaix, Walker, G. A. H., Walker, A. R., Irwin, A. W., Larson, A. M., Yang, 

D., Apps, K., & Laughlin, G. 2002, ApJ, 568, 352 S. L. S., & Richardson, D. C. 1995, Icarus, 116, 359 

Walker, G. A. H., Bohlender, D. A., Walker, A. R., Irwin, A. W., 
Yang, S. L. S., & Larson, A. 1992, ApJ, 396, L91 




Fig. 1. — The dependence of mean internal i) ra( j error on the stell^ $J<Fjbbtcd rotational velocity, vsini. The values of vsini and their 
sources are listed in Table 1. Vertical error bars indicate the standard deviation from the mean for all of our velocity measurements for a 
given star. 
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Fig. 3. — Comparison of our computed S index with published values from Duncan ct al. (1984). The straight line is the linear fit discussed 
in section 2.3. Error bars indicate the rms about the mean S value for each star. 
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Fig. 5. — Comparison between measured tr' v for our stars, and predicted cr' v based on empirical based on relations (see SBM) between a' v 
and P ro t (top left), vsini (top right) and R' HK (bottom left; separate fits for G and K stars used). Bottom right shows comparison between 
crj, (Hyadcs) and those predicted by the simple spot/plage model in SBM (smaller symbols have estimated P ro t)- In each case, F, G, and K 
stars are given by *, +, and diamonds, respectively, and the solid line indicates the observed a' v = a{,' prod . 
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Table 1 

Internal Errors versus Rotation for the Program Stars 



HD 


Other Name 


v sin i 
[km s -1 ] 


[m s" 1 


Reference 
] 


18632 


BD+07 459 


1.8±1.0 


3.52±0.72 


1 


26756 


vB 17 


6.4±3.0 


4.94±0.69 


2 


26736 


vB 15 


7.0±0.7 


5.13±0.99 


3 


27282 


vB 27 


6.6±3.0 


3.81±0.45 


2 


27406 


vB 31 


12. Oil. 


6.19±1.22 


3 


27859 


vB 52 


8.0± 0.7 


5.76±0.80 


3 


20899 


vB 64 


5.0±0.7 


3.81±0.46 


3 


28205 


vB 65 


9.0±3.0 


7.02±0.77 


4 


28237 


vB 66 


8.0±3.0 


6.80±0.90 


4 


28344 


vB 73 


7.0±0.7 


6.80±1.39 


3 



References. — (1) Fckel (1997);Wc adopted an error of 1.0 km 
s _1 . (2) Strassmeier et al. (2000); We adopted errors of 3.0 km 
s _1 . (3) Sodcrblom (1982) (4) Kraft (1965); Wc adopted errors 
of 3.0 km s _1 . 



Table 2 

Comparison of Observed and Predicted/Modeled a' v 



Type of 
Fit/Model 1 


Spec, type 
studied 


N 


K (prod -> -<) 

[dex] 


°Yms 

[dex] 


&rms (Lick) 2 

[dex] 


«hk(F) 


F 


11 


0.124 


0.444 


0.41 


«hk(F) 


F 3 


10 


0.210 


0.376 


0.41 


s hk( g ) 


F 3 


10 


-0.107 


0.213 




«hk(G+K) 


F 3 


10 


-0.213 


0.284 




r'(G+K) 


G 


24 


-0.045 


0.216 




R' HK (G+K) 


K 


34 


0.117 


0.348 




R' HK (G+K) 


K 3 


33 


0.146 


0.320 




fl^ K (G+K) 


G+K 


58 


0.050 


0.298 


0.25 


«hk(G+K) 


G+K 3 


57 


0.065 


0.278 


0.25 


s hk( g ) 


G 


24 


0.073 


0.224 




«hk(K) 


K 


34 


-0.040 


0.329 




-Rhk( k ) 


K 3 


33 


-0.011 


0.284 




v sin i(F) 


F 


3 


0.125 


0.263 


0.31 


t;sini(G+K) 


G 4 


7 


-0.016 


0.265 


0.22 


Prot(F) 


F 


1 


-0.008 




0.35 


P rot (G+K) 


G+K 


14 


0.021 


0.176 


0.22 


model 


F+G 4 


5 


0.001 


0.121 


0.14 


model 5 


F+G 4 


7 


0.013 


0.130 


0.16 



1 c.g., i?^ K (G+K) indicates that the cr / ( pred ) based on the empirical fit for G 
and K stars of a' v (Lick) vs. R' HK (from SBM) is compared to Hyades stars of 
the spectral type listed in column (2). Note: separate G star (oc #hk ) an< ^ 
K star (oc R^ 1 ' 19 ) fits for R' HK based on Lick data were not listed in SBM. 

2 from SBM 

3 with one outlier removed 

4 No Hyad K dwarfs in our sample have published v sin i measurements 
5 Assuming fs=0 for F stars 
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Table 3 
Program Information 



HD or HIP # 


Other Name 


# obs. 


B-V 




(7 <r'> 

, LLL 


r 




K(r 2 ) 


un 1/11 o n 
HD 14127 


BD+04 378 


12 


0.567 


2.84 


0.38 


0.55 


0.07 


0.66 


H1H lobUU 


r i r i i i-7 a c rr 
BD+17 455 


10 


0.704 


1.85 


0.31 


-0.08 


0.83 


2.41 


H1H looUo 


vB 153 


18 


0.855 


4.18 


0.49 


0.09 


0.73 


3.16 


HD 18b32 


BD+07 459 


9 


0.926 


4.00 


0.52 


0.00 


1.00 


2.33 


UT'i i nnnn 

hd iyyu2 


tz> n i on con 
HD + oz 582 


10 


0.732 


2.67 


0.40 


0.07 


0.85 


2.43 


MD ZU4oU 


vB 1 


13 


0.567 


2.07 


0.37 


0.62 


0.02 


0.25 


HD 20439 


vB 2 


13 


0.617 


2.88 


0.31 


0.30 


0.33 


1.79 


TITD 1 CCCQ 

H1.T lOODO 


nn i n7 /ino 
LSD + Ur 4yy 


9 


1. 130 


3.89 


0.38 


-0.60 


0.09 


0.62 


Hl.r 15 ( zU 


' ' ' 


8 


1 .431 


2.40 


0. 18 


0. 10 


0.82 


2. 14 


H1H lb529 


vB 4 


9 


0.844 


4.25 


0.36 


0.13 


0.73 


2.21 


tttt> i cnno 

hih ibyuo 


vB 5 


10 


0.917 


4.04 


0.38 


0.32 


0.36 


1.77 


HL) 23453 


nn i n rr n i o 

HD +25 bio 


8 


1.437 


2.43 


0.22 


-0.27 


0.51 


1.87 


H1H 177bb 


G7-15 


7 


1.340 


2.69 


0.29 


0.02 


0.96 


2.03 


TTTD 1 oni O 

H1H 18018 


vB 170 


11 


1. 160 


3.06 


0.36 


0.30 


0.38 


1.79 


TTn noCQCQ 
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